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We discuss the behaviour of the energy spectrum of particles in jets near the limit 
of small momenta of a few hundred MeV. In QCD parton cascades the soft gluons 
are coherently emitted from all faster partons in the jet and their production rate 
is predicted to scale. The observed charged and identified particle spectra follow 
this behaviour surprisingly well supporting of the hypothesis of the Local Parton 
Hadron Duality (LPHD) for this extreme limit. Further tests of this perturbative 
approach are discussed. 



1 Introduction 

The study of multiparticle production in hard collision processes can yield 
valuable information about the characteristic features of the partonic branching 
processes in QCD and the transition from the coloured partons to the colourless 
hadrons. The parton branching process is driven by the gluon bremsstrahlung 
which becomes singular in the limit of collinear or soft emission. Analytical cal- 
culations in the Double Logarithmic Approximation (DLA) take into account 
the contributions from these singularities. An improved accuracy is obtained 
in the Modified Leading Logarithmic Approximation (MLLA) which includes 
corrections of relative order y/al (for a review, seed). 

There is no commonly accepted theory of hadronization, but experimental 
data can be used to improve our knowledge about the partonic interactions 
in the limit o£ small momentum transfer. Here we consider the hypothesis 
of the LPHD □ which assumes that sufficiently inclusive observables at the 
hadronic level are directly given by the corresponding quantities computed for 
the parton cascade. A striking success of this approach is the prediction of the 
energy distribution of particles with an approximately Gaussian shape in the 
variable £ = log Ej e t /E for particles..with energy E in a jet of energy Ej et (the 
so-called "hump-backed plateau" ptatl. 

Here we focus on the soft end of the spectrum (E < 1 GeV). We recall that 
the gluons of long wave length are emitted coherently by the total colour cur- 
rent which is independent of the internal structure of the jet and is conserved 



1 



when the partons split. The soft radiation then depends essentially on the 
total colour charge of the initial partons. Due to the coherence of the gluon 
radiation it is not the softest partons but those with intermediate energies 
[E ~ -Ej e t~°' 4 ) which multiply most effectively in QCD cascades. Applying 
the LPHD hypothesis one expects that the hadron spectrum at low momen- 
tum p should be nearly independent of the jet energy i? je p&El. Quantitatively, 
we consider the invariant density Edn/cPp = dn/dyd 2 pT in the limit of van- 
ishing rapidity y and transverse momentum px, or, equivalently, for vanishing 
momentum p, i.e. 

„ dn 1 , „ dn 
I = lim £- r = -lim£- r 1 

y^Q :PT ^Q d A p 2 p^O d 6 p 

where the factor 1/2 takes into account that both hemisphere are added in the 
limit p — > 0. If the dual description of hadronic and partonic final states is 
really adequate down to very small momenta, the finite, energy independent 
limit of the invariant hadronic density Io is expected from the colour coherence 
argument. This expectation is supported by the experimental data for both 
charged and identified hadrons in the full energy range explored so far in e + e~ 
annihilation. Furthermore we discuss how the QCD picture explaining this 
result can be further tested. For more details we refer to the paperQ. 



2 Particle production in the soft limit 

We consider first the analytical predictions for the energy spectrum of partons 
near the soft limit. The asymptotic behaviour of the energy spectrum is ob- 
tained from the DLA in which energy conservation is neglected and only the 
leading singularities in the parton splitting functions are kept. The evolution 
equation of the single parton-inclusive energy distribution originating from a 
primary parton A is given bjtl: 

D p A (tY) = 6? A 6(0+£dt> £ ^ dy'^(y')DP(e,y')- (2) 

Here we have used the logarithmic variables £ = log(l/a;) = \og(Q/E) and 
Y = log(<2/<2o) with E the particle energy and Q the jet virtuality (Q = PQ for 
a jet of primary momentum P and opening angle 0); C a is the respective colour 
factor, Nc for A = g and Cf for A — q; 70 denotes the anomalous dimension 
of multiplicity and is related to the QCD running coupling by 7q = ANcol s J2-k 
or 7 g = /3 2 /log(p ± /A) with (3 2 = AN c /b, b = (11A C - 2n/)/3; A is the QCD- 
scale and Nc and n/ are the number of colours and of flavours respectively. 
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The shower evolution is cut off by Qq, such that the transverse momentum 

Px > Qq- n 

In case of fixed a, the exact DLA solution is knowm 



(3) 



where 7i is a modified Bessel function of first kind. This distribution in £ 
vanishes for £ — > and £ — > X. as required by eq. (ph. The exact MLLA 
correction is obtainedB (see alsola) by multiplying the DLA result (0) with the 
factor exp(-a 7 2(r - Q/W c ) where a = + |$r. 

In case of running a s the first terms of an iterative solution of the DLA 
equation (0) have been derived; the MLLA correction is again given by an 
exponential damping factor in analogy to the fixed a s caseQ. 

While the limiting behaviour of the partonic energy spectrum in the soft 
region follows from the general principle of colour coherence the detailed form 
of the observable hadronic spectrum is predicted uniquely from the LPHD 
hypothesis only for E ~ p 3> Qo, but not near the kinematical boundary 
because of the sensitivity of the spectrum to the cut-off procedure. In the 
usual application of the MLLA the partons are treated as massless with energy 
E = p > p± > Qq, so £ < Y. Experimental hadronic spectra are usually 
presented as function of momenta p or £ p = log(l/a;p) which is not limited 
from above. The same kinematic limit for partons and hadrons is obtained if 
the hadronic mass mh and the partonic p± cut-off Qo are taken the same. 

For the relation between parton and hadron distributions one may require 
that the invariant density Edn/d 3 p of hadrons approaches a constant limit for 
p — > as is observed experimentally. For the spectra which vanish linearly as 
in (J3J) this is achieved by relating the hadron and parton spectra asum 

Eh— -j = K h E p — = K h D y A {£ E , Y) (4) 

dph ap p 



with Eh — \fpl + <3o — E p > Qq, where Kh is a normalization parame- 
ter referring to a single jet. Then, indeed, for hadrons the invariant density 
Edn/d 3 p = K h D 9 A (£_, Y)j^{E\ - Q§) approaches the finite limit as in ([[]) 



3 



1 



[GeV- 



v/s = 91 GeV" 




,/s = 3 GeV 



MLLA 
DLA 



E[GeV] 

Figure 1: Invariant density Edn/dPp as a function of the particle energy E for Qq = 270 
MeV. Predictions of MLLA and DLA with fixed ct s (with 70 = 0.64) at cms energies of 
t/s = 3 GeV (lower two curves) and 91 GeV (upper two curves) with Dg computed using 
eq. $1). 



In the fixed a s limit /3 2 /A is replaced by 70 and Iq ~ 1/Qo- With prescription 
(||) the moments of the full energy spectrum D(£, Y) are well described by the 
MLLA formulae at Qo — 270 MeV in a wide energy range&D. 

The relation (0) is not unique, however. We found that the alternative 



prescription based on phase space arguments, dn/dt; 
e.g.0): 



(p/E) 3 D(Z,Y) 



d 3 Ph 



AttE 2 



Y) 



(G) 



works well for charged pions for all energies E in the LEP region and for 
charged particles at low energies E if D 9 A is the MLLA "limiting spectrum" 
with Qo = A = 138 MeV. The low cut-off mass is plausible in this region, 
which is dominated by pions. 

In order to illustrate the above exact analytical results we compare in 
Fig. 1 the computations from DLA and MLLA (in case of fixed a s ) for low 
particle energies using the relation (|]) between parton and hadron spectra. 
The single particle invariant density approaches an energy independent value 
in the soft limit £ — > Y (E — > Qo)- This originates from the soft gluon emission 
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Figure 2: a) Invariant density Edn/d?p of charged particles in e+e - annihilation as a 
function of the particle energy E = \J "p 1 + Qq at Qo = 270 MeV at various cms energies 
with Dg computed using eq. (^) (K^ = 0.45); b) the same as in a), but at Qo = 138 MeV 
with Dg computed from the Limiting Spectrum using eq. (H) (K^ = 1.125). 

contribution of order a s (the term of order 7q in the expansion (^)) which is 
determined by the total colour charge of the primary partons due to the colour 
coherence. In this limit the MLLA converges towards the DLA as the energy 
conservation constraints are unimportant and the parton splitting functions 
are only probed for very small fractional momenta. The approximate results 
with running a s show the same features!! 

3 Discussion of experimental results 

Fig. 2a shows the charged particle invariant density, Edn/d 3 p, as a function 
of the particle energy E at different cms energies ranging from 3 GeV up to 
LEP-1.5 cms energy (133 GeVjUun comparison with the MLLA predictions for 
running a s in the approximate forntJ. In this figure, eq. (^) is used to relate 
partor^and hadron spectra at the effective mass Qo = 270 MeV which was 
found U to provide a good description of the moments of the energy spectra 
over a large cms energy interval. In Fig. 2b the same data are compared 
with the MLLA limiting spectrum (A = 0), but using eq. (0) with the mass 
value Qo = 138 MeV. Whereas this mass value is plausible for small energies 
E, the fit does not work well for the full region of energies E or the lower 
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cms energies. Note that, while these fits behave differently very close to the 
boundary E ~ Qq, they clearly demonstrate the scaling behaviour expected 
from colour coherence. Furthermore they reproduce the considerable energy 
dependence of the initial slope which can be related to the running of a s in 
the MLLA damping factor. 

It is remarkable that the data from all cms energies tend to converge in 
the soft limit; we find I ~ 6-8 GeV~ 2 (using Q = 270 MeV) and I ~ 4-6 
GeV~ 2 (using Qo = 138 MeV). Inspecting the soft limit more closely, the LEP 
data seem to tend to a limiting value larger by about 20% as compared to the 
lower energy data. This may be well due to the overall systematic effect in the 
relative normalization of the different experiments. Alternatively, there could 
be a contribution from incoherent sources, for example, from weak kaon and 
heavy quark decays. I-. 

The good scaling behaviour of the soft limit Iq is also found in case of 
identified particles n, K , and p. Remarkably, the scaling continues down to 
the lowest energy where data are available, at y/s=1.6 GeVlla. 

4 Further tests of the parton hadron duality picture 

It is suggestive to relate the observed limiting behaviour of particle spectra 
for vanishing momenta to the expected behaviour of soft gluon emission off 
the primary partons. The particle rate Iq in this limit (see eq. (^)) cannot 
be directly predicted, as it depends on the normalization factor and the cut- 
off parameter Qq. However, the remarkable scaling behaviour of Iq in e + e~ 
annihilation provides one with a standard scale for the comparison with other 
processes. It would be very interesting to study whether the soft particle 
production is universal, i.e., a purely hadronic quantity or whether the intensity 
Jo indeed depends on the colour topology of the primary active partons in the 
collisions process. This could provide one with a direct support of the dual 
description of soft particle production in terms of the QCD bremsstrahlung. 

In what follows we consider some possibilities to test this hypothesis further 
and discuss what one may expect from the comparison of the e + e~ reaction 
with other collision processes. 

4-1 gg Final State 

A direct test of the relevance of QCD coherence for the soft production limit is 
the comparison of the qq with the gg colour singlet final state. In the gg final 
state the same line of argument applies as above for the qq final state, but the 
colour charge factor is increased by Nc/Cf = 9/4 in eq. (||) implying roughly 
a doubling of the soft radiation intensity Iq- An approximate realization of 



6 



the colour octet antenna is possible in e + e — * qqg with the gluon recoiling 
against a quasi-collinear qq painj. 

4-2 S-Jet Events in e + e~ Annihilation 

To be more quantitative, we consider the soft radiation in 3-jet qqg events 
of arbitrary jet orientation into a cone perpendicular to the production plane 
and compare it to the radiation into the same cone in a 2-jet qq event again 
perpendicular to the primary qq directions. The analysis of these configurations 
avoids the integration over the kr > Qo boundaries along the jets. We restrict 
ourselves to calculations of the soft gluon bremsstrahlung of order a s . From 
the experience above with the two-jet events this contribution dominates in the 
soft limit whereas higher order contributions take over with increasing particle 
energy. 

Let us consider first the soft radiation, into arbitrary direction ft from a qq 
antenna pointing in directions fti and nJi-3: 



dN qq = ^-dQ H7 ^W^(H) , W™(H)=2C F (fj) (7) 



with (ij) — ciij /(diCij), a,ij = (1 — niftj) and ctj = (1 — ririi). Such an antenna is 
realized, for example, in a qqj event, and can be obtained by the appropriate 
Lorentz boost from the qq rest frame. The soft gluon radiation in a qqg event 
is given as in (^) but with the angular factor 

W™(n) = N C [(1+) + (M - j^(=P)] (8) 

where (+,—,1) refer to (q,q,g). 

For the radiation perpendicular to the primary partons (ij) — aij = 1 — 
cos6,y, with the relative angles 0^ between the primary partons i and j. For 
2-jet events of arbitrary orientation one obtains in this case 

Wf(<d+-) = 2C F (l - cos9 + _) . (9) 

Correspondingly the ratio R± of the soft particle yield in 3-jet events to that 
of 2-jet events in their own rest frame (W^ 9 (7r) = 4Cf) is given by 

dNf 9 N c 1 

R± = -^r = - cos i+ - cos e i- - jvtC 1 ~ cos ©+-)]■ ( 10 ) 

This yields, in particular, the limiting cases R±_ — 1 for soft or collinear primary 
gluon emission (0i+ = 1 — ©i-, ©h = 7r) and the effective gg limit R± = 
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Nc/Np for the parallel qq (0+_ = 0) configuration, as expected. Eq. (10) 
also predicts R± for all intermediate cases, in particular, for Mercedes-type 
events one finds R± = l.59. In this case no jet identification is needed for the 
above measurements. We also note that the large angle gluon radiation is 
independent of the mass of the quark (for O ^> mq/Ej et ). 



4-3 Hard Processes with Primary Hadrons and Photons 

We consider here processes with dominantly 2-jet final states. These include 
semihard processes which are initiated by partonic 2-body scatterings. In case 
of quark exchange (in particular, .deep inelastic scattering processes 7yp, jvl 
at large Q 2 in the Breit frame Ii3't£l) the two outgoing jets originate from the 
colour triplet charges and Iq should be as in e + e~ annihilation; in the case of 
gluon exchange Iq should be about twice as large (Nc/Cp) as expected for a 
gg jet system. In 2-jet events the limiting density Iq may be most conveniently 
determined as a limit of dn/dyd 2 pT at y ~ for px — > 0. In general the deep 
inelastic processes may have contributions from both quark exchange and gluon 
exchange which are expected to dominate in different frames. Final states with 
several well separated jets can be treated in analogy to the e + e~ — > 3-jet events 
discussed before. 

The frame dependence can be studied conveniently in these exchange pro- 
cesses by considering Io(y) — dyd i pT | PT _> as a f unc ti° n °f rapidity y, mea- 
sured, say, in the cms frame. The Breit frame is reached by a Lorentz transfor- 
mation in the direction of the incoming photon, and corresponds to a different 
rapidity, y Breit, in the cms. One would then expect for the processes with the 
virtual photons at large Q 2 in general a "quark plateau" of Io(y) in the cur- 
rent region of length Ay ss log(Q/m) near y Breit and a transition to a "gluon 
plateau" in the complementary region of length Ay « \og(W 2 / Qm) where m 
is a typical particle mass. The existence of a plateau corresponds to the energy 
independence of Iq as seen in e + e~ annihilation. So, if gluon exchange occurs 
with sufficient hardness in the process, Iq should develop a step like behaviour 
but never become larger than twice the Iq value in e + e~ annihilation. The 
height of the "plateau" is a direct indicator of the underlying exchange process 
(quark or gluon like). 



4-4 Soft Collisions (Minimum Bias Events) 

These processes (with initial hadrons or real photons) are not so well under- 
stood theoretically as the hard ones but it might be plausible to extrapolate 
the gluon exchange process towards small pjO. Experimental data at ISR 
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Figure 3: Invariant density Edn/d?p in hadronic collisions at central rapidity as a function 
of pt at different cms energies. 

energies, however, do not follow the expectation of a doubling of Iq in com- 
parison to e + e~ annihilation, rather the intensities are found to be similaii3. 
On the other hand, as shown,, in Fig. 3, Iq roughly doubles when goingJrom 
- 20 GeV at SppS/ISRliaka to ^ = 900 GeV at the SppS CollideM If 
additional incoherent sources like weak decays can be excluded, such behaviour 
could indicate the growing importance of one-gluon exchange expected from 
the perturbative picture. Then a saturation at Iq H /Iq e ~ 2 is expected and 
no additional increase for a semihard pr trigger. A rise of Iq h could also result 
from incoherent multiple collisions of partons (e.g. E3). 

4-5 Production of W pairs in e + e~ annihilation 

If each of the W's decays into 2 jets the limiting soft particle yield will be 
twice the yield in e + e~ — > qq , provided the interconnection phenomena^ are 
neglected. This tests directly the hypothesis of independent emission. 

5 Conclusions 

The analytical perturbative approach to multiparticle production, based on the 
MLLA and the LPHD, has proven to be successful in the description of various 
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inclusive characteristics of jets. It is of importance to investigate the limitations 
of this picture, in particular in the soft region where non-perturbative effects 
are expected to occur. 

It is remarkable that the intensity of the soft hadron production follows 
to a good approximation a scaling law in a range of two orders of magnitude 
in the cms energy of e + e~ annihilation (1.6-140 GeV). Such a scaling law is 
derived analytically for the soft gluons in the jet and follows directly from 
the coherence of the soft gluon radiation from all emitters. It appears that 
the production of hadrons which is known to proceed through many resonance 
channels nevertheless can be simulated in the average through a parton cascade 
down to a small scale of a few 100 MeV as suggested by LPHD. The scaling 
behaviour down to the very low cms energy of 1.6 GeV can only be explained 
within a parton model description if the cut-off scale Qo is well below 1 GeV. 

It will be interesting to investigate the suggested scaling property (1) and 
its possible violation further in the same experiment to avoid systematic effects; 
this seems to be feasible at LEP and HERA for the simplest processes. In 
particular, the effect of weak decays on scaling violations should be clarified. 
The sensitivity of the soft particle production to the effective colour charge of 
the primary emitters can be tested through the transverse production rates 
in multijet events. The soft radiation in e + e~ annihilation can be used as a 
standard scale in the comparison of various processes. In this way the soft 
particle production can be a sensitive probe of the underlying partonic process 
and the contributions from the incoherent sources. 
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